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Long-range exchange and correlation effects, responsible for the failure of currently used approximate density 
functionals in describing van der Waals forces, are taken into account explicitly after a separation of the electron- 
electron interaction in the Hamiltonian into short- and long-range components. We propose a "range-separated 
hybrid" functional based on a local density approximation for the short-range exchange-correlation energy, 
combined with a long-range exact exchange energy. Long-range correlation effects are added by a second-order 
perturbational treatment. The resulting scheme is general and is particularly well-adapted to describe van der 
Waals complexes, like rare gas dimers. 
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I. INTRODUCTION 

Van der Waals (dispersion) interactions are universal attrac- 
tive forces due to long-range correlation of electrons between 
weakly- or non-overlapping electron groups 0]. They play an 
important role in the cohesive energy of practically all kinds 
of materials: intermolecular complexes, extended systems, 
like molecular crystals, liquids or biological macromolecules. 
Although, in principle, density functional theory (DFT) 0] 
within the Kohn-Sham (KS) scheme |3[ is able to provide 
the exact ground state energy of an electronic system, present 
approximate density functionals are inappropriate to describe 
long-range electron correlation and consequently fail for van 
der Waals interactions, manifested by their incapacity of re- 
producing the correct R~ e asymptotic behaviour of the inter- 
molecular potential 0] ■ 

Several propositions have been published recently to add 
the missing long-range correlation contribution or to use 
asymptotically correct correlation energy expressions in DFT. 
Most of these methods require a partitioning of the system 
into interacting parts and are valid only for large separa- 
tions 0,0, 01- Seamless dispersion energy functionals 0], 
that are valid for the whole range of possible intermolecular 
separations have also been proposed 1 9] . A general problem in 
schemes that use an additive correction to standard function- 
als is the double counting of a part of the correlation effects 
that are already present in the original functional. 

Moreover, it is not enough to add missing correlation ef- 
fects to traditional density functionals. Many of the present 
approximate functionals, like the local density approximation 
(LDA) which is well-known for its notorious overbinding ten- 
dency, and also many popular generalized gradient approx- 
imations (GGA), already predict a more-or-less pronounced 
bound state for simple van der Waals complexes, like rare gas 



dimers 1 10]. As it has been pointed out by Harris twenty years 
ago lllll . this behaviour is related to the erroneous distance 
dependence of approximate exchange functionals. In effect, 
in self-interaction corrected calculations the minimum on the 
potential curve disappears J5]]. Therefore, in order to describe 
correctly both the minimum and the asymptotic region of van 
der Waals potential energy surfaces it is mandatory to remove 
the unphysical bonding by appropriately correcting the ex- 
change functional. 

Here, we propose a scheme based on a long-range/short- 
range decomposition of the electron interaction which meets 
the above requirements and remedies the description of van 
der Waals forces in the framework of a first-principles ap- 
proach, which takes into account simultaneously long-range 
correlation and exchange effects, avoids double counting and 
is size-extensive. 

Our scheme is based on the hypothesis that for the descrip- 
tion of van der Waals (London) dispersion forces one should 
improve the representation of long-range electron interaction 
(exchange and correlation) effects. At a first level of approx- 
imation, we treat the long-range exchange energy explicitly 
while maintaining a density functional approximation for the 
short-range exchange-correlation energy. This step defines a 
"range-separated hybrid" (RSH) scheme, which is corrected 
in a second step for the long-range correlation effects by a 
second-order perturbation theory, leading to size extensive 
M0ller-Plesset (MP2)-like correction. This method will be re- 
ferred to as RSH+MP2. 

The idea of a long-range/short-range decomposition of the 
electron interaction is not new (see, e.g., Refs. [3, IB Q3, 
In the context of DFT, this approach has been 
used to construct multi-determinantal extensions of the KS 
scheme (HEEHJElEIHHi. A density functional 
scheme with correct asymptotic behaviour has been proposed 
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along these lines very recently by Baer and Neuhauser ll25ll 
and the correct 1/r asymptotic behaviour of the long-range 
exact exchange has been also exploited in time-dependent 
DFT calculations of polarizabilities constituting the 

major motivation of the recent development of "Coulomb at- 
tenuated" hybrid functionals |3- Heyd, Scuseria and Ernz- 
erhof applied an inverse range-separation in order to get rid of 
the convergence problems of the exact exchange in solid-state 
calculations I29l 13011 . Their HSE03 functional is a generaliza- 
tion of the PBEO hybrid functional 1I3TI1 where the long-range 
portion of the exact exchange is replaced by the long-range 
component of the PBE exchange functional 13211 . 

In the context of the calculation of van der Waals ener- 
gies, the idea of separating the electron interaction operator to 
short- and long-range components has already been explored 
by the work of Kohn, Meier and Makarov, who applied the 
adiabatic connection - fluctuation-dissipation approach for 
long-range electron interactions l33ll . leading to an asymptot- 
ically correct exp ression of the dispersion forces. It has also 
been shown 13411 that the artificial minimum of the rare gas 
dimer potential curves can be removed by an exact treatment 
of the long-range exchange. 

The second order perturbational treatment of the full 
Coulomb interaction has already been used by several authors 
for the van der Waals problem 1I5 H35H36I1 . and it was shown 
that the resulting asymptotic potential has the qualitatively 
correct 1/R 6 form. As shown very recently, quantitatively 
reliable asymptotic form of the potential energy curve can be 
expected from adiabatic connection - fluctuation-dissipation 
theory calculations 13711 . 

The general theoretical framework is outlined in SectionlTTl 
describing the RSH scheme and the second-order perturba- 
tional treatment of long-range correlation effects. As de- 
scribed in Section [HI] our approach has been tested on rare 
gas dimers. These systems are typical van der Waals com- 
plexes, where the attractive interactions are exclusively due 
to London dispersion forces. They constitute a stringent test 
of the method, since the potential curves have very shallow 
minima of the order of about 100 /jM. 

Unless otherwise stated, atomic units is assumed through- 
out this work. 



n. THEORY 

A. Multi-determinantal extension of the Kohn-Sham scheme 

We first recall the principle of the multi-determinantal ex- 
tension of the KS scheme based on a long-range/short-range 
decomposition (see, e.g., Ref . uM and references therein). 

The starting point is the decomposition the Coulomb 
electron-electron interaction w ee (r) = 1/r as 

w ee (r)=w*f(r)+w^(r), (1) 

where w^(r) = erf (/ir)/r is a long-range interaction and 
w ee^ ( r ) i s me complement short-range interaction. This de- 
composition is controlled by a single parameter fi. For [i = 0, 



the long-range interaction vanishes, w^ =0 (r) = 0, and the 
short-range interaction reduces to the Coulomb interaction, 
w ee M=0 (r) = w ee (r). Symmetrically, for /1 — > 00, the 
short-range interaction vanishes, w s ^^°° (r) = 0, and the 
long-range interaction reduces to the Coulomb interaction, 
u, ee M ^°°( r ) = w ee{r)- Physically, 1/fi represents the dis- 
tance at which the separation is made. 

The Coulombic universal density functional F[n] = 
min*^ n (*|f + Weel*} 03], where T is the kinetic en- 
ergy operator, W ee = (1/2) JJdr 1 dr 2 w ee (r 1 2)n2(r 1 ,r 2 ) is 
the Coulomb electron-electron interaction operator expressed 
with the pair-density operator 7i2(ri, r^), is then decomposed 
as 

F[n\ = F^[n]+E^ c [n\, (2) 

where F lr ^[n] = miri*^„(*|T + W^f\%>) is a long- 
range universal density functional associated to the interac- 
tion operator = (1/2) jjdr 1 dT2w\f{r 1 2)fi2{r l , r 2 ), 
and £hxcM = ^'"N + E£f[n] is by definition the cor- 
responding complement short-range energy functional, com- 
posed by a trivial short-range Hartree contribution ' M [n] = 
(1/2) jj dYidY2W^(ri2)n(Yi)n(Y2) and an unknown short- 
range exchange-correlation contribution E^[ri\. At /i = 0, 
the long-range functional reduces to the usual KS kinetic en- 
ergy functional, F k,A1=0 [n] = T s [n], and the short-range func- 
tional to the usual Hartree-exchange-correlation functional, 
EftJ/ c ~°[n] = -Ehxc M ■ In the limit fi — > 00, the long- 
range functional reduces to the Coulombic universal func- 
tional, F 11 "'^^ 00 [n] = F[n), and the short-range functional 
vanishes, E^ ( T'° [n] = 0. 

The exact ground-state energy of a A^-electron sys- 
tem in an external nuclei-electron potential w„ e (r), E = 

Yiiu\ n ^N$yF[n] + Jdrw ne (r)ri(r)| where the search is over 

all A^-representable densities, can be re-expressed using the 
long-range/short-range decomposition of F[n] 

E= mint F b *[n] + E%£.[n] + /dru„ e (r)n(r)} 

n—*-N I J J 

= minf (*|f + W^/l*} 

*^JV I 

+ + «M}. (3) 

where the last search is carried out over all A^-electron nor- 
malized (multi-determinantal) wave functions , 3>. In Eq. Q, 

(r) is the density coming from the wave function i.e. 

(r) = (\f r |n(r)|\I r ) where h(r) is the density operator. 
The minimizing wave function "f^ in Eq. (0 is given by the 
corresponding Euler-Lagrange equation 

(f + W*f + Vne + KfcM) I*") = £"1*"), (4) 

where V ne = J drw„ e (r)n(r), V^cM = J rfr %^( r )"( r ) 
with the short-range Hartree-exchange-correlation potential 

v Hxc( r ) = ^Hic[ n ]/^ n ( r )' anc l ^ i s me L a g ran g e multi- 
plier associated to the constraint of the normalization of the 
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wave function. Eq. (|4) defines a long-range interacting effec- 
tive Hamiltonian H» = f +W**+V ne +V££ [n^] that must 
be solved iteratively for its multi-determinantal ground-state 
wave function \I r/i which gives, in principle, the exact phys- 
ical ground-state density n(r) = n*f»(r) = (4' AI |ri(r)|5' AI }, 
independently of fi. Finally, the exact ground-state energy ex- 
pression is thus 

E = (*"|f + + V ne \V) + £££[n*,]. (5) 

This exact formalism enables to combine a long-range wave 
function calculation with a short-range density functional. In 
the special case of \i = 0, the KS scheme is recovered, while 
the limit fj, — ► oo corresponds to the usual wave function for- 
mulation of the electronic problem. 

A short-range LDA 01 and other beyond-LDA 01 
approximations have been constructed to successfully de- 
scribe the functional [n] . In previous applications of the 
method, the long-range part of the calculation has been han- 
dled by configuration interaction |2lll or multi-configurational 
self-consistent field (MCSCF) 1 23] methods. We propose in 
this work to use instead perturbation theory. 

B. Range-separated hybrid 

At a first level of approximation, we introduce the RSH 
scheme by restricting the search in Eq. l|3} to iV-electron nor- 
malized one-determinant wave functions $ 

£ M ,rsh = min (,$>f + 

+ /drv„ e (r)n$(r) +£^[n$]j. 

(6) 

The associated minimizing one-determinant wave function 
$ /i satisfies the Euler-Lagrange equation 

(f + v ne + + v&*,]) m = wm, 

(7) 

where V h '^hf[ < I ) ] is a long-range potential operator appearing 
due to the restriction to one-determinant wave functions 
as in Hartree-Fock (HF) theory, and £q is the Lagrange 
multiplier associated to the normalization constraint. As 
usual, V^. M HF [$] is the sum of a Hartree contribution, 

v£&m n dv 1 dv 2 w b e f(r 12 mh(v 1 )mh(v 2 ), 

and a non-local exchange contribution, V^'hfI^] = 
-(1/2) II dr 1 dr 2 w^{r 12 )(^\ni^2 7 r 1 )\^)n 1 {ri,r 2 ), 
where rii(ri,r2) is the first-order density matrix operator. 
Eq. defines the RSH non-interacting effective Hamiltonian 
&g = f + V ne + V^ HF [^] + V*f c M that must be 
solved iteratively for its one-determinant ground-state wave 
function $ M . Of course, does not give the exact physical 
density: ^n. 

The RSH energy expression is finally 

E^ Rm = <#"[T + V ne \^) + + E^ c [n^\, 

(8) 



where E^ n¥ [<5>] = (<5>\W { ^\<5>) is the HF-like long-range 
Hartree-exchange energy. Eq. defines a single-parameter 
hybrid scheme combining a long-range HF calculation with 
a short-range density functional. The case fj, = still corre- 
sponds to the KS scheme while the method reduces now to a 
standard HF calculation in the limit \i — ► oo. 

We note that an equivalent to the RSH scheme has been 
investigated recently by Pedersen and Jensen 12311 as a special 
case of the combination of a long-range MCSCF calculation 
with a short-range density functional. 

C. Long-range correlation corrections by perturbation theory 

We develop now a long-range perturbation theory, using the 
RSH determinant $ p as the reference. To do so, we introduce 
the following energy expression with a formal coupling con- 
stant A 

E^ x = min{(*|f + Vne + V^ x %[^] +AW lr '^|*) 

+*££M}. (9) 

where the search is carried out over all iV-electron normalized 
(multi-determinantal) wave functions and W 1 '^ is the long- 
range fluctuation potential operator 

W"'" = - Vh^hp^]. (10) 

The minimizing wave function tyv> x j n Eq. ||9} is given by 
the Euler-Lagrange equation 

(f + v ne + v^ HF m + aw 1 ^ + v^cK^l) l*"' A ) 

(11) 

where E^' is the Lagrange multiplier associated to the nor- 
malization constraint. For A = 1, the physical energy is re- 
covered, E = E^' 1 , in principle independently of n, and 
Eq. <|TT) reduces to Eq. (g): 1^ A=1 = £^ A=1 = £/*. 
For A = 0, Eq. (II Q reduces to the RSH effective Schrodinger 
equation of Eq. 0: <I"'< A=0 = £M. A =o = £ p 

We expand E» x in powers of A, E^ x = J^kLo E^X k , 
and apply the general results of the non-linear Rayleigh- 
Schrodinger perturbation theory lEH l42l l43ll outlined in the 
Appendix. It is easy to verify that the sum of zeroth- and first- 
order energy contributions gives back the RSH total energy 

The second-order correction can be written as 

(13) 

where Rq is the reduced resolvent 
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in terms of the excited eigenfunctions <f>^ and eigenvalues £q l 

of the RSH effective Hamiltonian Hjf, and Gq is a short-range 
screening operator 



2 |y drdr'n(r)|^)/^[n*,](r,r')(^|n(r'), 

(15) 

with the short-range Hartree-exchange-correlation kernel 
£W(r,r') = 6 2 E^ c [n]/5n(r)Sn(r>). 

Let insert the spectral resolution of Eq. dl4i in Eq. Jl 3I >. 
Since W 1 ' is a two-electron operator only singly and doubly 
excited determinants, and $!"-„,, where i, j refer to oc- 
cupied spin-orbitals and a, b to virtual spin-orbitals of can 
a priori contribute to _E^>( 2 ) . Actually, singly excited deter- 
minants gives vanishing matrix elements, ($f_^JW lr |<I> M } = 
0, since it can be easily verified that (^f^jW^I^) = 
($f^al^Hx M HF[ $/ 1l* M )' as in standard HF theory. Conse- 
quently, the product RqGq in Eq. (II 31 involves vanishing 
matrix elements, ($' i |n(r)|$^_ >ab ) = 0, i.e. the non-linear 
terms are zero with the present choice of the perturbation op- 
erator W lr . The second-order energy correction is thus 



£M,(2) 



E 

a<b 

E 

a<b 



m^ a b\w^m\ 2 



CM 
c 



(16) 



where is a spin-orbital of $ M and is its associated eigen- 



value, 



are the two-electron integrals asso- 



ciated to the long-range interaction w^(r 12 ), and we recall 
that the indexes i, j refer to occupied spin-orbitals and a, b to 
virtual spin-orbitals. Eq. (I16> is fully analogous to the conven- 
tional MP2 energy correction. The total RSH+MP2 energy is 

^/x,RSH+MP2 _ £jn,RSH _|_ J7jp,(2) 

^From a practical point of view, once the RSH orbitals and 
one-electron eigenvalues are available, any standard MP2 im- 
plementation can be used, provided that the long-range elec- 
tron repulsion integrals corresponding to the RSH orbitals are 
plugged in. Due to the long-range nature of these integrals one 
can take advantage of efficient modern algorithms, like the lo- 
cal MP2 H, multipolar integral approximations, which have 
particularly favorable convergence pr oper ties for long-range 
part of the split Coulomb interaction 14511 . or the resolution 
of identity approach j46Tl . It means that in appropriate im- 
plementations the extra cost of the MP2 corrections can be 
made negligible for large systems with respect to the resolu- 
tion of the self-consistent RSH equations, similar to a usual 
KS calculations with a hybrid functional. Solid state applica- 
tions for semi-conductors can also be envisaged on Wannier 
orbital-based MP2 implementations l47ll . 



System 


r/™ (a u ^ 


£™ (//FT) 


C« (a u ) 


C xp (a u ) 


He 2 


5.62 


34.87 


1.534 


1.461 


Ne 2 


5.84 


134.18 


6.860 


6.282 


Ar 2 


7.10 


454.50 


73.19 


63.75 


Kr 2 


7.58 


639.42 


153. 1 


129.6 



TABLE I: Absolute parameters of the reference potential curves de- 
termined from Ref. |53 The C 6 ' coefficients were obtained from a 
logarithmic fit in the same conditions as explained for the calculated 
potentials. 



III. RESULTS AND CONCLUSIONS 

The above described RSH+MP2 approach has been applied 
to rare gas dimers, using a LDA-based short-range exchange- 
correlation functional with a range-separation parameter of 
/i = 0.5. This latter value corresponds to the smallest mean av- 
erage error of the atomization energies calculated by the RSH 
scheme for the G2-1 set (a subset of 55 molecules of the G3 
set 1 48, |0]) of small molecules l5(ill . This value is in agree- 
ment with the intuitive picture predicting that should be 
close to the physical dimensions of a valence electron pair. 
The interaction energies were calculated with a modified ver- 
sion of the MOLPRO package 15 ill . The basis set superpo- 
sition error (BSSE) has been removed by the counterpoise 
method. 

The results are presented as reduced potentials, U*(r*) = 
U(r* -d^/sm, where the reduced variables U* = U/e m 
and r* = r/d m , are defined with respect to the equilibrium 
distance d m and the well-depth e m of accurate "experimen- 
tal" potential curves l52ll (cf. Table The calculated po- 
tentials are characterized by the hard core radius, a* defined 
by U*(cr*) = (experimentally a* « 0.89), the reduced well 
depth, U* n , and the equilibrium distance, r* n , (experimentally, 
by construction, U* n = — 1 and r* n — 1). The minimum region 
is also characterized by the harmonic vibrational frequencies, 
co, related to the second derivative of the potential at the min- 
imum. 

The long-range behaviour of the potential energy curves 
can be appreciated from the C6 coefficients. Experimen- 
tal C6 coefficients are usually obtained from optical data 
(dipole oscillator strength distributions) 115311 and characterize 
the purely dipolar contribution to the long-range interaction 
energy. Since we had no access to such a decomposition of the 
interaction energy, we have determined an effective Cq coef- 
ficient by a logarithmic fit of the interaction energies between 
30 and 60 Bohrs. This quantity, which includes higher order 
multipolar effects too, is presented in the form of a reduced 
variable, Cg = Cq/Cq 1 . Here C| l has been obtained from an 
analogous fit to the points of the reference potential reported 
in Table|I] For the sake of comparison, the experimental C 6 xp 
(purely dipolar) values are also reported. 

The RSH and RSH+MP2 potential curves, as well as the 
HF, the standard MP2 and the coupled-cluster CCSD(T) ones, 
calculated with the aug-cc-pVTZ basis set are represented for 
the four dimers in Figure ^ and compared to the experimen- 
tal curves. Note that the reduced representations of the ex- 
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He dimer Ne dimer 




0.8 1.0 1.2 1.4 1.6 0.8 1.0 1.2 1.4 1.6 



Ar dimer Kr dimer 




0.8 1.0 1.2 1.4 1.6 0.8 1.0 1.2 1.4 1.6 



FIG. 1: Reduced HF (dotted repulsive), RSH (dashed repulsive), 
MP2 (long dashed), CCSD(T) (dashed), RSH+MP2 (dotted) and 
Tang-Toennis reference (full) potential curves for He2, Ne2, Ar2 and 
Kr2 dimers. 



perimental potentials of different rare gas dimers are prac- 
tically indistinguishable. The calculated RSH potentials are 
always repulsive, like the HF ones. The RSH+MP2 poten- 
tials are slightly too repulsive at short interatomic distances, 
as reflected by the values of the hard core radii, systematically 
higher than the experiment (around 0.89). The RSH+MP2 
and CCSD(T) curves are almost the same for Ne2 with a well 
depth of around U* n = 0.6, while the RSH+MP2 minima of 
the Ar 2 and Kr2 systems are even better (U^ > 0.9) than the 
CCSD(T) ones (U* n ps 0.7). The position of the minimum 
is predicted within 1-4% in the RSH+MP2 approximation. 
The 6-8% deviation found for the He 2 RSH+MP2 minimum 
can be explained by an exaggerated repulsion, reflected by the 
highest a* found in this case. In comparison with the usual 
MP2 potential curves, the RSH+MP2 follow similar trends, 
being systematically more stabilizing and closer to the exper- 
imental curve. 

The main quantitative features of the RSH+MP2 poten- 
tials obtained by the aug-cc-pVTZ and aug-cc-pV5Z basis 
sets 154115 51 15«. 15711 are summarized in Table Inland compared 
to the results of standard MP2 and CCSD(T) supermolecule 
calculations with the same basis sets. The basis set has a 
non-negligible effect on the calculated parameters of the po- 
tential curves, which converge systematically towards the ex- 
perimental values for all the properties. For He 2 the double 
augmented d-aug-cc-pV5Z basis set results are also included, 
representing a further improvement of the well depth, but hav- 
ing practically no effect on the equilibrium distance. 

The basis set superposition error of the equilibrium dis- 
tances and of the interaction energies are reported in Tablelllll 
as the difference in the parameters of the BSSE-contaminated 



and BSSE-free reduced potential energy curves. The BSSE 
corrections on the bond lengths and on the interaction ener- 
gies are always negative, i.e. the BSSE-contaminated dis- 
tances are too short and the energies are too low. In some 
case, like the Ne 2 dimer with aug-cc-pVTZ basis, the bind- 
ing energy correction may attain 55 or 67% of the well depth 
at the MP2 and CCSD(T) level of approximation. The cor- 
responding RSH+MP2 BSSE effect is considerably smaller, 
but it is still 34%. The BSSE effect on the bond lengths are 
much less spectacular, but still more pronounced in the MP2 
and CCSD(T) methods than in the RSH+MP2 approach. As 
a general trend we can conclude that the RSH+MP2 has usu- 
ally less than the half of the MP2 or CCSD(T) basis set su- 
perposition errors. This is a considerable advantage for an 
efficient and reliable exploration of potential energy surfaces, 
especially when the lack of well-defined subsystems make im- 
possible to perform a counterpoise correction. 

Effective Q coefficients obtained from the RSH+MP2 ap- 
proach agree with the experiment within 5% for He 2 and 
Ne 2 , and are overestimated by 15-20% for Ar 2 and Kr 2 . It 
means that the asymptotic behaviour of the RSH+MP2 poten- 
tial curves is reasonable. We recall that the exact Cq coeffi- 
cient is given by the Casimir-Polder relation Ql 

C(j = — / G?wai(z<x>)a 2 (io;) (17) 
n Jo 

where a\ (iui) and a 2 (iui) are the exact dynamical polarizabil- 
ities of the monomers. It is known that the asymptotic form 
of the MP2 energy expression corresponds to an uncoupled 
HF-type, non-interacting approximation of the monomer po- 
larizabilities. This means that MP2 calculations do not repro- 
duce the exact Cg coefficients: usually they tend to overes- 
timate them. For instance, in the case of the benzene dimer, 
this overestimation in the complete basis limit may reach a 
factor of 2; for less polarizable systems the situation is less 
critical. An analogous behaviour is expected for RSH+MP2. 
Note however, that in this case one-electron excitations are 
obtained from the self-consistent RSH one-electron states, 
which include, in addition to the long-range exact exchange, 
short-range exchange-correlation effects too. A more reli- 
able approximation can be developed on the basis of the adi- 
abatic connection - fluctuation-dissipation approach l33l l37ll 
which would ensure, in principle, the exact asymptotic limit 
of the potential energy curves. The development of a range- 
separated version of this method is under progress. 

In conclusion, the RSH+MP2 approach provides an ef- 
ficient DFT-based description of weak intermolecular com- 
plexes bound by dispersion forces. Even in its simplest, LDA- 
based implementation, it represents a huge improvement over 
KS calculations, which lead to unreliable potential curves 
in the minimum region with a qualitatively wrong asymp- 
totic behaviour. Range-separated extensions of other density 
functionals, like the gradient-corrected PBE functional, are in 
progress. By removing systematic errors of currently used ap- 
proximate DFT functionals and introducing corrections which 
grasp the essential physics of van der Waals interactions, the 
RSH+MP2 approach extends the applicability of density func- 
tional calculations to weak intermolecular forces. Further 
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Method MP2 CCSD(T) RSH+MP2 



System 


' m 


U* 




a 


Ce 


r m 


U* 


u/uo m 


a 


Ce 


r m 


u; n 




a 


Ce 


He 2 
AVTZ 
AV5Z 
d-AV5Z 


1.000 
1.052 
1.036 
1.032 


-1.000 
-0.516 
-0.594 
-0.629 


1.000 
1.058 
1.032 
1.032 


0.894 
0.941 
0.926 
0.923 


1.000 
0.754 
0.760 
0.763 


1.000 
1.023 
1.007 
1.003 


-1.000 
-0.777 
-0.896 
-0.946 


1.000 
1.047 
1.019 
1.020 


0.894 
0.912 
0.896 
0.893 


1.000 
0.966 
0.980 
0.958 


1.000 
1.080 
1.078 
1.077 


-1.000 
-0.553 
-0.593 
-0.613 


1.000 
1.141 
1.135 
1.135 


0.894 
0.961 
0.957 
0.955 


1.000 
1.008 
1.028 
1.038 


Ne 2 

AVTZ 

AV5Z 


1.000 
1.073 
1.043 


-1.000 
-0.435 
-0.588 


1.000 
1.004 
0.977 


0.896 
0.960 
0.936 


1.000 
0.766 
0.816 


1.000 
1.041 
1.009 


-1.000 
-0.609 
-0.877 


1.000 
0.983 
0.950 


0.896 
0.931 
0.904 


1.000 
0.914 
0.990 


1.000 
1.040 
1.036 


-1.000 
-0.605 
-0.751 


1.000 
0.965 
0.965 


0.896 
0.928 
0.923 


1.000 
0.950 
1.036 


Ar 2 

AVTZ 

AV5Z 


1.000 
1.023 
0.998 


-1.000 
-0.850 
-1.062 


1.000 
1.033 
0.996 


0.897 
0.913 
0.891 


1.000 
1.095 
1.136 


1.000 
1.037 
1.011 


-1.000 
-0.715 
-0.910 


1.000 
1.051 
1.013 


0.897 
0.927 
0.903 


1.000 
0.958 
1.175 


1.000 
1.012 
1.007 


-1.000 
-0.948 
-1.040 


1.000 
1.013 
1.003 


0.897 
0.903 
0.896 


1.000 
1.154 
1.215 


Kr 2 

AVTZ 

AV5Z 


1.000 
1.029 
1.002 


-1.000 
-0.840 
-1.080 


1.000 
1.049 
1.021 


0.896 
0.918 
0.893 


1.000 
1.124 
1.153 


1.000 
1.048 
1.016 


-1.000 
-0.677 
-0.898 


1.000 
1.069 
1.038 


0.896 
0.936 
0.908 


1.000 
0.960 
0.980 


1.000 
1.016 
1.007 


-1.000 
-0.919 
-1.023 


1.000 
1.011 
1.008 


0.896 
0.905 
0.897 


1.000 
1.136 
1.154 



TABLE II: Reduced parameters of the calculated MP2, CCSD(T) and RSH+MP2 (/i =0.5) potential energy curves obtained by the aug-cc- 
pVTZ (AVTZ), aug-cc-pV5Z (AV5Z) and d-aug-cc-pV5Z (d-AV5Z) basis sets. Reduced experimental parameters are listed in the first line for 
each dimer. Absolute reference values are given in TableQ] 



Method 


MP2 


CCSD(T) 


RSH+MP2 


System 


' 772, 


EC 


* m 


m> 


' m 




He 2 














AVTZ 


-0.007 


-0.125 


-0.008 


-0.121 


-0.000 


-0.063 


AV5Z 


-0.004 


-0.048 


-0.002 


-0.037 


-0.001 


-0.015 


d-AV5Z 


-0.008 


-0.160 


-0.007 


-0.113 


-0.001 


-0.031 


Ne 2 














AVTZ 


-0.049 


-0.547 


-0.035 


-0.674 


-0.031 


-0.335 


AV5Z 


-0.011 


-0.150 


-0.006 


-0.148 


-0.001 


-0.025 


Ar 2 














AVTZ 


-0.020 


-0.263 


-0.023 


-0.239 


-0.007 


-0.101 


AV5Z 


-0.005 


-0.138 


-0.004 


-0.103 


-0.002 


-0.022 


Kr 2 














AVTZ 


-0.013 


-0.191 


-0.017 


-0.174 


-0.007 


-0.126 


AV5Z 


-0.003 


-0.073 


-0.002 


-0.049 


-0.002 


-0.039 



from 7i* (r) = (\&|n(r)|\E f ), where n(r) is the density op- 
erator. In Eq. JAIL A is a formal coupling constant; we are 
ultimately interested in the case A = 1. The minimizing wave 
function ty x satisfies the Euler-Lagrange equation 



(ff (0) + XW + S1 A ) 



= £ x \* x ) 



(A2) 



where the eigenvalue £ A comes from the normalization condi- 
tion and Cl x is a potential operator coming from the variation 
of F[n], non linear in A, 



dv 



8F[r 



Sn(r] 



(A3) 



TABLE III: BSSE correction for the reduced parameters r* m and t/JJ 



tests should decide whether this method is generally applica- 
ble to the important domains of the physisorption, or cohesion 
in molecular crystals and in layered semi-conductors. 



APPENDIX A: NONLINEAR RAYLEIGH-SCHRODINGER 
PERTURBATION THEORY 

Let us consider the following general total energy expres- 
sion, involving a Hamiltonian H^°>, a perturbation operator 
W and a density functional F [n] , 



{(*|2T(°) +\W\*)+F[rii S ,}}, (Al) 



E = min 

*->iV 



where the search is carried out over all jV-electron normalized 
wave functions "f, (^l^) = 1, and n* is the density coming 



where n x is the density coming from ty x , n\(r) = 
(* A |n(r)|* A ). 

Starting from the reference A = 0, we develop a perturba- 
tion theory in A. We introduce the intermediate normalized 
wave function ty x 



l* A > 



(A4) 



and expand 4' A , n A , fl x and £ A in powers of A: ^ A = 

IXo ^ (fc)Afe ' nX = EZo n(k)xk > ^ A = Er= ^ {k)xk and 

£ A = J2T=o £ (fc) A fe . The coefficients are obtained from 
the expansion of ^ A through 



n A (r) = 



(* A |n(r)|^ A ) 
(* A |§ A ) 



(A5) 



and the coefficients fl^' are found from the expansion of n A , 
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after expanding tl x around rS°\ 

Q A = / dr „ 1 ; ; J n(r) 



Sn(r) 

dvdv' f F }f?\ An\r>)h(r) + --- . 
on{r)dn{r') 

(A6) 

where An A = n A — n^ ) . The zeroth-order equation is 

|*<°>> = £(°)|^°)), (A7) 
and of course = v]/ A = . For the general order k > 1, 

(AW + - £W) |*W> + VKI*^- 1 ') 

k k 

+ ^0«|^( fe -*)) £(*)!*(*-*)>. (A8) 

i=l i=0 

The corresponding eigenvalue correction of order fc is 

k 

£( fc ) = (|r(0)|^|f (fe-1)) + ^(f (0)|QW|^r(fc-i)) ) (A9) 

i=l 

containing, besides the usual first term, a "non-linearity" term 
as well. Introducing the reduced resolvent, Rq, 



\¥ I 0) )(¥° ) \ 

J" £<?> - £(°) 



(A10) 



where ¥^ > * > and £^ are the excited eigenfunctions and eigen- 
values of H(°\ the wave function correction of order k writes 

\¥ k} ) = -R w\¥ k - 1] ) - R Ci {k} \¥ a] ) 



fc-i 



-R 



(^ ( '° ~ £(i) ) |* (fc-<) )- (AH) 



The total energy can be re-expressed in terms of the eigen- 
value £ A and the "double counting correction" D x 

E x = £ A + D X , 



(A12) 



where 



D x = F[n x ] - [ dr^^-n x (r). (A13) 
J x n{r) 

We expand E x and D x in powers of A: E x = J^kLo E {k) X k 
and D x = £^ .D (fc) A fe . The coefficients are found 
from the expansion of n x , after expanding D x around rS°\ 

D x = F[n<°>] + [ dr^^An A (r) 
J on(r) 

Iff.., 5*F[n(°)} A _ A .,. A . 



I drdr 6n(r)6n(r>) 



^- An (r')An (r) + 



dr 



SF[n 



(0)1 



V(r) 



<5n(r) 

1 . An (r )n (r) — ■ ■ ■ 
dn(r)on(r') 



drdr' 



The zeroth-order total energy is simply 

#(o) = £ (o) + F[n ( 0)l _ / dr ^V!H n (o )(r); (A15) 



<5n(r) 

The general correction of order k > 1 writes 

= (^r(O) |^|^r(fc-l)) + A( fc ) (A16) 

where A*^ is 

fe 

A (fc) = E(* (0) |fi (<) l* (fc_i) > +D^ k l (A17) 
i=l 

At first order, it can be verified that the nonlinearity term of 
the eigenvalue and the double counting correction cancel each 
other, i.e. A^ 1 ) = 0, and we obtain the conventional first-order 
energy correction 



(A18) 



At second order, the situation is analogous, i.e. A' 2 ' = 0, 
and again the conventional form of the energy correction is 
retrieved 



£( 2 ) = (¥ ^\W\¥^) 



(A19) 



The nonlinearity effects are "hidden" in the first-order wave 
function correction, which can be obtained from the self- 
consistent equation: 

ifC 1 )) = -R W\¥ Q) ) - £ Q (1) \¥ Q) ) (A20) 

Since the first-order potential operator is, for real wave func- 
tions, 



dn(r)dn(r) 

(A21) 

Eq. (IA20> can be re-expressed as 

= -R W\¥ 0) ) - i? Go|* (1) ), (A22) 



(A14) 



where 

Go = 2 // d rdr'h(v)\¥^) /f^ (¥^\n(r'). 
J J dn(r)dn(r) 

(A23) 

The final expression of the second-order energy correction can 
be written as the series 

#(2) = -(¥^\W(l + RoG )~ 1 R W\¥ '>) 
= -(¥V\WRoW\¥ 0) ) 
+ (¥^\WRoG R W\¥^} 

. (A24) 

Further details and higher-order expressions will be given in a 
forthcoming publication. 
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